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Abstract 

We consider the exact solution of a model of correlated particles, which is pre- 
sented as a system of interacting XX and Fateev-Zamolodchikov chains. This 
model can also be considered as a generalization of the multiband anisotropic 
t — J model in the case we restrict the site occupations to at most two elec- 
trons. The exact solution is obtained for the eigenvalues and eigenvectors 

using the Bethe-ansatz method. 

PACS numbers: 75.10.Lp, 74.20-z, 71.28+d 
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It is well known that in the limit of a strong Coulomb repulsion U the traditional Hubbard 
model reduces to the so-called t — J model [1]. The strong on-site correlations limit the 
site occupations to at most one electron. States with double occupation on a given site are 
energetically unfavorable and can be projected out from the Hilbert space. In one dimension 
this model is integrable at the supersymmetrical point [2-5]. The exact solution can be 
obtained also for the extension of the supersymmetric spin-1/2 t — J model (i.e. N = 2) to 
the case of arbitrary number TV of bands, having a SU(N) symmetry [6,7], as well as for 
their anisotropic version [8,9]. In some sense all these models can be considered as a system 
of interacting XX and XXZ chains with size occupations limited to at most one electron per 
site [2]. 

In this paper we consider a generalization of the multiband t — J model for the case 
where the on-site correlations are enough strong to limit the site occupations to at most 
two electrons. The exact integrable model we will present can be considered as a system of 
interacting XX and spin-1 Fateev-Zamolodchikov chains [10,11]. 

In order to present our model let us consider initially a general model with N distinct 
type of bands (a = 1, . . . , N) split into two disjoint groups Ki and N 2 . Taking into account 
the constraint of maximum double occupancy in a given site we can have at most one 
electron in bands a G Ki and at most two electrons in bands a G K 2 . The states of a given 
site are denoted by \ct, (3 > (a < (3, a, (3 — 0, 1, . . . , N), where in the case a, (3 — 1,2, ... N 
these are the bands where the electrons are located while a = or (3 = denotes the 
absence of electrons. As a conequence of this notation the state \a,a > is forbidden if 
a G Ni and it is allowed if a G K 2 or a = 0. The most general Hamiltonian with nearest- 
neighbours interactions, in a lattice with L sites and periodic boundary condition, that 
conserves separately the electrons in each band can be written as 

L 

N 

TT Tl/CCl, 02,03,^4 p Q l' Q 2|oi,a2 JTi«3,a4| 03,04 /-,x 

[a],[a']=0 
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where [a] = [a[, a' 2 , a' 3 , a' 4 ] is a permutation of [cui, «2, 0:3, 04] which keeps the order 

ttl < «2, tt'l < «3 < «4, « 3 < 04, Ctj = 0, 1, 2, N. 

The matrix £ ,a '^l a ''' 3 ' has all the elements zero except the element in the a, (5 line and ol \(3' 
column, with unit value, i. e. 

E^'* 13 ' = \a[3><a'{3'\. (2) 
We consider Hermitean Hamiltonians, 

W a ,/3 ' 7 ' s — W a ' l3 '' y ' 5 ' (1) 

in the absence of external fields, 

wSsgjg = <o ;o°;« = o, o < a < n, (4) 

and satisfying the chirality property 

W^j, (rj) = W^^i-v) , (5) 

where 77 is the asymmetry parameter. We choose the energy scale such that all single-particle 
hopping couplings have a unit value, 

WSg$$ = t = l 0<a<N. (6) 

Let us consider initially the case of a single particle on the otherwise empty chain. As 
a consequence of translational invariance of the Hamiltonian (1) the eigenf unctions are the 
plane waves of wavenumber k, 

*i = E/(^«iKT l|0 ' l >---> >> 

XI 

f(x 1 ,a 1 ) = e ikx , k=—l (Z = 0,1,..,L-1), (7) 

with energy 

E = -2 cos A;. (8) 
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In (7) |0, . . . , > is the reference state with no particles. 

The eigenfunctions of the Hamiltonian (1) in the general case, where we have n particles, 

^n = EE «i5 ^ a 2 ; ...; x B , a n )£°f ^'Xf »IW ..i£f»l >. . . . , >, (9) 
{Q} M 

are calculated by using a generalized nested-Bethe ansatz. In (9) the first summation is 
over all the permutations Q = [Qi, Q 2 , ■ ■ ■ , Q n ] of the integers 1, 2, . . . , N, and for a given 
permutation [Q] the second sum is restricted to the set 1 < xq 1 < xq 2 < • • • < xq u . We 
seek the amplitudes / in each of these regions in the form of a superposition of plane waves 
with wavenumbers kj,j — 1, . . . ,n. If we have only single occupation (xQ i 7^ XQ i+1 ,i = 
1, 2, . . . , n — 1) we write the ansatz 

n 

f(x u a Ql ;x 2 , a Q2 ; x n , a Qn ) = ^ ^pf.!."p n Q " II ex P( ik Pj x Q 3 )> ( 10 ) 

p 3=1 

where the sum is over all permutations P — [Pi, . . . , P n ] of the integers 1, 2, . . . , n. In the 
case we have a pair at the position xq l = xq 1+1 , the ansatz is modified to 



f(xi,a Ql ;x 2 , a Q2 ; x n , a Qn ) = ^ ^ZpX+i-^p^ II «^i}k Pj x Qj ), (11) 

p i=i 

where the bar at the /th and (/ + l)th positions of the superscript indicates the pair location. 
The general case with many isolated particles and pairs follows from (10) and (11). The 
ansatz (10)-(11) corresponding to the configurations where \xQ i+1 —XQ i | > 1 satisfies H\^f >= 
E\fy >, providing the energy and momentum are given by 

n n 

E = -2 ]T cos kj\ P = J2 k J- ( 12 ) 
3=1 3=1 

In order to obtain the two-particle scattering matrix we consider the general amplitudes 
(10) or (11) in the case where there are only two particles on two neighboring sites. In this 
case our problem is exactly equivalent to that of a model with only two species (bands). 
The coefficients A^ ' P ™ Qn arising from the different permutation Q are connected with each 
other by the elements of the two-particle S-matrix 

N 

A:tp,.. = E s^(k Pl ,k P M::P 2Pl ... (is) 

5,7=1 



As a necessary condition to ensure integrability, the two-particle scattering matrix has to 
satisfy the Yang-Baxter equation [12 - 13]. There are two integrable models, described in 
terms of two species of particles, whose S-matrix has a factorizable form. The first of these 
models is the standard Hubbard model [13], while the second one is the correlated hopping 
model [14 - 16]. We may try to use as the fundamental building block of a general S-matrix 
of our model (1) the two-particle scattering matrix of these models. Here we consider the 
second model and choose the S-matrix as in the correlated hopping model [16] []. In this 
case we have the following restrictions on the parameters of the Hamiltonian (1) 

<5S? = 0, 0<a<(3 

W3gS? = eo 0<a<{3 

W to°'° = ICS? = v / iT^ ) < a < (3 

Wojafija = 2ei cosh rj, a G K 2 

W0S8g = eo 0<a<(3 } (14) 



where eo = ±1 and e% = ±1. The non- vanishing elements of the S-matrix are 



ToPfj. j„ \ _ sin(Ai-A 2 ) 



sin(Ai - A 2 - irf) 

laPfu u ^ _ Zl sinh 1 

sin(Ai - A 2 - irj) 



Sfa(k 1} k 2 ) = ^t^t^ — Yexp[zsign(a - /3)(Ai - A 2 )] 



SSfe.^)-- 8 " 11 "^ 1 "^"''' 1 (15) 
sm(Ai — A 2 — irj) 



where 



sin( Xj + ir]) 

e lk] = -e — VT 2 rr- (16) 



1 Unfortunately previous attempts to solve this problem by using the Hubbard S-matrix was 
unsuccessful [17] since the Bethe ansatz does not work in the sector where we have three or more 
particles. The interactions necessary to prevent more than double ocupancy per site spoils the 
exact integrability. This is not the case if we use as the fundamental building block the S-matrix 
of the correlated hopping model. 



To complete the proof of the Bethe ansatz (10-11) we must consider the eigenvalue equa- 
tions in the case where there are three and four particles on the two neighboring sites j and 
j + 1. This gives us a complicated system of equations for those parameters of the Hamil- 
tonian (1) involving three and four particles. We have treated this system analytically and 
checked our results numerically. We found the following solution for the diagonal elements 
of W 



W^ = e (l-e^), 0< 7 < a < f3, 
WSSS? = 0. 0<a< 7 </3, 



e , a G Ki 



e e 2r? , aGH 2 



W£& 4 = eo(l - 2e 2 "), 0<a</^<7<5, 

KlSS = -e e 2 \ 0<a< 1 <[3<6, 

w Zt$ = e o(l - 2cosh2?7), < a < 7 < 5 < (3, 
eo(l-e^), 13 



-2e e 2 ", (3 G K 2 

.eoll-e 2 ''), 7^1 
W2$Z={ > 0<a</?< 7 , 



-2eo cosh2?7, 7 G H2 
e (l-e 2 "), aeNi 



— 2e cosh2r7, a G K 2 
2e , a^eNi 



— 2e cosh 277, ct or /3 G K 2 



KSSl = - e o(2e 2r? + e" 2 "), a G K 2 , < a < [3, 



= + 2e^), (3 G K 2 , < a < (3, 

WZZ = \w%%£ = -2e cosh 277, a G K 2 . (17) 



For the nondiagonal elements of W we have 
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<?#7 = 1 0<a</3< 7 , 

<ft 7 7 = Wg^f = e~\ < a < < 7 , 

, 2 cosh 77, fleNi 
W'Zo Z = < , < a < /3, 



0, /3eK 2 
2 cosh?7, a G Ki 



' 0, a G K 2 

<'£o /I = = eiv / TT^, < a < /3, 

K3$$ = <fe 7 = < a < < 7 , 

= - £ o, < a < (3 < 7 , 



wSZSg = wSSSS = -^i vTT^, < a < /?, 

WSJ?? = «4J = "^o, 0<a</?< 7 <5, 

a,/?,7,<5 rv a,-y,l3,8 

= W^W^ S = -e e\ < a < /? < 7 < 5, 

= W^W^ = -e e 1 VTT^, < a < /3 < 7, 

{— 2e cosh 77, /9 G Ni 
0, /? G N 2 

WKSf = <SS? = -60, < a < 0, 

= "2eo cosh rj < a < (3, 
Wl$$ = *> <*,P*>y,6. (18) 

The periodic boundary condition for the system on the finite lattice, with size L, gives 
us the Bethe-ansatz equations. In order to obtain these equations we must diagonalize the 
transfer matrix of a related inhomogeneous vertex model with non-intersecting strings [18]. 
The Bethe-ansatz equations are written in terms of the momenta of the electrons kj and 
additional rapidities 
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Bm(Aj + irj) L = £n _ mi _! -p, sin(A j - \ r + iejrj) ^ sin (A, - A£> - fay) 
[ sin(A,-*7?) 1 jUi si^Aj - Ajv - ir/) ^ sin (A j - A« + 

?} sin(AW - AiT } + t eg+1 7y) = nv _ 1+f?lff _ 1 ^+^-1 "ff sin(AW - Ag^ + |^) 
ii sin(A^ - A<? - fe^) CT a+1 ii sin^ - A^ X) - je^iy) 

„ ^ sin(AW - A^ 1 ) + i gg+1 „) 
' ii sin(AL CT) - A^ +1) - | £(T+1 r,) ' 

a = 1,2, ...,N- 1, m Q = n, m N = 0, A^ 0) = A a , 

where rij = mj_i — is the number of particles of species j, e a = +1 if a G Ki and e CT = — 1 
for a G H 2 . The energy of the system is given in terms of the Bethe-ansatz roots Xj 

E = -2gcos/, = 2 eo g[cosh2,- cosh2>; _J s2Aj ]. (20) 

It is interesting to remark that the particular isotropic cases 7 = 0, of (19 - 20), where N = 2 
and £1 = — e 2 = 1, as well as for arbitrary N, but £j = l(i — 1, . . . , N) have been obtained 
in refs. [19] and [20], respectively. 

In conclusion, we have studied a new integrable model which can be presented as a 
system of interacting XX and Fateev-Zamolodchikov chains. The Bethe-ansatz equations 
are derived by means of the coordinate Bethe-ansatz approach. A desirable continuation 
of the present work is the investigation of the thermodynamic equilibrium properties of the 
model based on the Bethe-ansatz solutions (19 - 20). 

This work was supported in part by Conselho Nacional de Desenvolvimento Cientifico e 
Tecnologico - CNPq - Brazil, by FINEP - Brazil and by the Russian Foundation of Funda- 
mental Investigation (Grant 99-02017646). 
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